UNCLASSIFIED 


_ AD  NUMBER _ 

AD867847 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative/Operational  Use;  FEB  1970.  Other 
requests  shall  be  referred  to  Chief  of  Naval 
Operations ,  OP-96 ,  Washington,  DC. 


_ AUTHORITY 

CNO  ltr  23  Oct  1970 


THIS  PAGE  IS  UNCLASSIFIED 


tfO  r 


NO  FORM 


LIEHABY 

TECHNICAL  REPORT  SECTION 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CAL1TCRHIA  93940 


AVAL  WARFARE  ANALYSIS  GROUP 


A  MODEL  OF  CARRIER-SUBMARINE 
INTERACTIONS 

By  J.V.  Hall,  Cdr.,USN 
J.K.  Tyson,  CNA 
J.S.  Finucane,  CNA 

CNA  Research  Contribution  No.  137 


V 


y 


This  Research  Contribution  does  not  necessarily 
represent  the  opinion  of  the  Department  of  the  Navy 


1401  Wilson  B  ou  lav  ard  Arlington.  Virginia  22209 


CONTRACT  N00014— 68— A  —  0091 


This  document  is  subject  to  special  export 
controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be 
made  only  with  prior  approval  of  the  Office 
of  the  Chief  of  Naval  Operations  (Op-96). 


* 


f 


U01  i>1 546 

pZscn^Amo°r:\ HaW S  Scant  J  IP  ALEXANDRIA  VA  N«AL  WARFARE  ANALYSIS  GROUP 
Report  Date:  February  01  1970 
Media:  39  Page(s) 

Distribution  Code:  01  -  APPROVED  FOR  PUBLIC  RFl  faqP 
Report  Classification:  Unclassified  LEASE 

Source  Code:  252950 
From  the  collection: 

Technical  Reports 


liO  FGRN 

Center 

for 


1401  Wilson  Boulevard  Arlington.  Virginia  22209  703/524-9400 


Analyses 


an  affiliate  of  the 
University  of  Rochester 


(NWG) 69-70 
15  April  1970 


From:  Director,  Naval  Warfare  Analysis  Group 

To:  Distribution  List 

Subj :  Center  for  Naval  Analyses  Research  Contribution  No.  137; 

forwarding  of 

Enel:  (1)  CNA  RC  137,  "A  Model  of  Carrier-Submarine  Interactions" 

by  CDR  J.  V.  Hall,  USN ,  J.  K.  Tyson,  and  J.  S.  Finucane 
Unclassified  dated  February  1970 


1.  Enclosure  (1)  is  forwarded  as  a  matter  of  possible  interest. 

2.  This  Research  Contribution  describes  a  parametric  model  of 
engagements  between  submarines  and  a  single  aircraft  carrier 
with  ASW  support.  The  measures  of  effectiveness  are  the  prob¬ 
ability  that  the  carrier  can  remain  on  station  for  a  specified 
length  of  time  and  the  expected  number  of  submarines  killed. 

The  model  was  developed  as  an  aid  in  the  study  of  how  various 
tactics  and  forces  can  improve  the  efficiency  of  carrier  operations. 

3.  Research  Contributions  are  distributed  for  their  potential 
value  in  other  studies  or  analyses.  They  have  not  been  examined 
in  detail  and  do  not  necessarily  represent  the  opinion  of  the 
Department  of  the  Navy. 


JAMES  K.  TYSON 
Acting  Director 


Distribution  List: 
Reverse  page 


(NWG) 69-70 
15  April  1970 


Subj :  Center  for  Naval  Analyses  Research  Contribution 

No.  137;  forwarding  of 


DISTRIBUTION  LIST 

OASD (SA) 

DIR,  WSEG  (2) 
ADMIN,  DDC  (20) 

CNO:  Op-09B9 

090 
095 
96 
31 

COMFIRSTFLT 
COMSECONDFLT 
COMSIXTHFLT 
COMSEVENTHFLT 
COMASWFORLANT 
COMASWFORPAC 
COMSUBPAC 
COMSUBLANT 
COMCARDIV  1 
COMCARDIV  2 
COMCARDIV  3 
COMCARDIV  4 
COMCARDIV  5 
COMCARDIV  6 
COMCARDIV  7 
COMCARDIV  9 
COMCARDIV  14 
COMCARDIV  16 
COMCARDIV  20 
COMASWGRU  3 
COMASWGRU  5 
COMSUBDEVGRU  TWO 


Op-32 

322 

05W 

06C 

71 


USNA  ANNA 

NAVPGSCOL 

NAVWARCOL 

NAVSHIPRANDCEN 

NAVAL  ASW  DATA  CENTER 

PM-4,  ASW  Systems  Project 


Other : 

APL/JHU 

IDA 

ORL 

RAND 

SRI 

SDC 

TRW  Systems  Group 

Cornell  Aeronautical  Laboratory 


CNA  RESEARCH  CONTRIBUTION  NO.  137 


NAVAL  WARFARE  ANALYSIS  GROUP 

CENTER  FOR  NAVAL  ANALYSES 


A  MODEL  OF  CARRIER -SUBMARINE 
INTERACTIONS 


By  J.  V.  Hall,  Cdr. ,  USN 
J.  K.  Tyson,  CNA 
J.  S.  Finucane,  CNA 


February  1970 


Work  conducted  under  contract  N00014-68-A-0091 


This  document  is  subject  to  special  export 
controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be 
made  only  with  prior  approval  of  the  Office 
of  the  Chief  of  Naval  Operations  (Op-96). 


Enclosure  (1)  to 

(NWG)69-70 

Dated  15  April  1970 


NAY  WAG  ASW  MODELS  AND  FORCES  STUDY 


K.H.  Bohlin,  Project  Director 

J.V.  Hall,  Cdr.,  USN,  Project  Officer 


ABSTRACT 


This  paper  develops  analytic  expressions  for  the  probability  that  a  single 
aircraft  carrier,  opposed  by  submarines,  can  remain  on  station  for  a  specified 
length  of  time.'  The  threat  from  both  torpedo  and  cruise -missile  submarines  is 
considered,  but  the  threat  from  aircraft  is  not.  Expressions  for  expected  sub¬ 
marine  losses  are  derived.  Expressions  are  also  developed  to  show  how  the 
probability  of  remaining  on  station  improves  as  the  carrier's  resistance  to  damage 
increases.  A  wide  variety  of  ASW  force  compositions  and  tactics  can  be  rep¬ 
resented  in  the  parameters  of  the  model.  Appendix  C  presents  an  alternative 
way  to  formulate  the  problem  and  derive  the  quantities  mentioned  above.  The  two 
methods  produce  results  that  are  in  excellent  agreement.  The  intention  of  the 
authors  is  to  provide  a  tool  that  can  be  used  to  improve  carrier  effectiveness 
through  the  study  of  tactics  and  force  interactions. 
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A  MODEL  OF  CARRIER-SUBMARINE  INTERACTIONS 


INTRODUCTION 

When  submarines  oppose  a  carrier  that  is  operating  in  the  same  general 
area  for  a  period  of  time,  the  events  that  take  place  are  very  complex  and 
cannot  be  exactly  described  in  mathematical  terms.  This  paper  defines  and 
models  a  simple  process  that  is  similar  in  some  respects  to  the  above  events. 
Suitably  interpreted,  the  model  can  be  used  to  study  the  interactions  of  a 
compact  or  dispersed  carrier  formation  and  torpedo  or  cruise-missile  sub¬ 
marines. 

The  following  rules  govern  the  modeled  process.  (1)  There  is  a  constant 
rate  of  encounter  1  (probability  of  encounter  per  unit  time)  between  the  carrier 
and  each  surviving  submarine.  (2)  At  each  encounter  between  the  carrier  and 
a  submarine  there  is  a  fixed  probability  that  the  carrier  will  be  successfully 
attacked  2  and  a  fixed  probability  the  submarine  will  be  killed.  (3)  As  long  as 
the  carrier  remains  on  station,  the  rate  of  attrition  to  each  surviving  submarine 
(probability  of  kill  per  surviving  submarine  per  unit  time)  is  constant.  When  the 
carrier  retires  all  action  stops  and  all  rates  become  zero.  (4)  The  submarines 
are  identical  and  their  operations  are  mutually  independent. 

The  following  notation  is  used  in  the  model: 

a  =  Rate  at  which  surviving  submarines  are  attrited  when  there 
are  no  successful  attacks  on  the  carrier  ("submarines  killed 
per  surviving  submarine  per  unit  of  time") 

b  =  Rate  of  encounter  between  each  surviving  submarine  and  the 
carrier  ("encounters  per  surviving  submarine  per  unit  time") 

h  =  Pr  (carrier  is  successfully  attacked,  given  an  encounter) 

a  =  Rate  at  which  surviving  submarines  are  attrited  in  the  absence 

of  encounters  with  the  carrier. 


Definitions  of  "encounter"  for  various  situations  are  given  in  the  "Discussion" 
section. 


A  "successful  attack"  is  an  attack  resulting  in  one  or  more  hits  on  the  carrier. 
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S  =  Unconditional  rate  at  which  surviving  submarines  are  attrited 
("submarines  killed  per  surviving  submarine  per  unit  time") 

N  =  Number  of  successful  attacks  the  carrier  can  absorb  before 
being  forced  to  retire 

x  =  Initial  number  of  attacker's  submarines 
o 

r  =  Rate  at  which  additional  submarines  arrive  in  the  area 

K(t)  =  Expected  number  of  submarines  killed  up  to  time  t 

R(t)  =  Ratio  of  expected  number  of  submarines  killed  to  the 

probability  that  the  carrier  has  been  forced  to  retire  up  to 
time  t 

P  (t)  =  Pr  (carrier  has  been  successfully  attacked  exactly  n  times 
n  by  time  t,  n=0,  1,  2,  N) 

tt  (t)  =  Pr  (carrier  has  been  successfully  attacked  n  or  fewer  times 
n  by  time  t,  n=0,  1,  2,  . . .  ,  N). 

The  first  8  quantities  listed  are  inputs  to  the  model;  the  last  4  are  outputs. 

SUMMARY  OF  EQUATIONS  DERIVED 


The  equations  summarized  in  this  section  are  derived  in  appendix  A. 


The  probability  that  the  carrier  has  not  been  successfully  attacked  by 
time  t  is 


PQ  (t)  =  exp  Jrt  +  (xq  -  r/a)  (l-e"at) 


(1) 


This  function  has  the  general  shape  shown  in  figure  1.  A  special  case  occurs 
when  r=0,  meaning  that  no  additional  submarines  arrive  in  the  area  after  the 
action  begins.  Then 


PQ(t)  =  exp 


bhx 


-at\ 
(1-e  ) 


(2) 


and  the  function  (shown  in  figure  2)  has  a  limiting  value  given  by 


lim  P  (t) 
o 

t  -»  co 


exp 


bhx 


(3) 


-2- 


Po  «> 


FIGURE  1 


P0  W 


o 


FIGURE  2 


-3- 


This  is  the  probability  that,  in  the  special  case  where  r=0,  the  carrier  will 
not  be  successfully  attacked  before  all  of  the  submarines,  x  in  number,  are 
killed.  0 

For  a  carrier  that  can  absorb  N  successful  attacks  before  being  forced 
to  retire,  the  probability  of  undergoing  exactly  n  successful  attacks  by  time 
t  is,  for  n  ^  N  , 


pn(t>  =  p0(t> 


[-In  P  (t)]  n 


n  ! 


(4) 


The  probability  that  such  a  carrier  has  undergone  n  or  fewer  successful 
attacks  is,  for  n  <N  , 


TT_(t)  =  E  Pk(0 
n  k=0  K 


(5) 


If  the  carrier  can  absorb  N  successful  attacks  before  being  forced  to  retire, 
its  probability  of  being  on  station  at  any  time  t  is  irN(t)  .  This  probability  is 

shown  in  figure  3,  plotted  against  1  -  P  (t)  ;  that  is,  figure  3  shows  the 

probability  that  the  carrier  is  still  on  station  versus  the  probability  that  the 
carrier  has  been  successfully  attacked  at  least  once.  The  feature  of  figure  3 
is  that  it  shows  the  relative  advantage  enjoyed  by  a  carrier  that  can  absorb 
some  successful  attacks. 


7TN  (t) 


FIGURE  3 
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by 


The  expected  number  of  submarines  that  are  killed  up  to  time  t  is  given 


K(t) 


_S_ 

bh 


N 

I  [l-rrk(t)]. 
k=0 


(6) 


The  "exchange  ratio"  is  given  by 


R(t) 


K(t) 

1  -  ttn  (t)  • 


(7) 


If  the  carrier  cannot  absorb  any  successful  attacks  and  must  retire  the  first 
time  it  is  hit,  then  N=0,  and  the  exchange  ratio  is  independent  of  the  variable 


R(t)  = 


_S_ 
bh  • 


(8) 


For  N  2:  1,  it  is  possible  to  calculate  the  expected  number  of  submarines  the 
enemy  would  have  to  expend  to  force  the  carrier  to  retire,  if  he  were  willing 
to  commit  submarines  for  as  long  as  necessary: 

lim  R(t)  =  (N+l)  -A.  .  (9) 

t  -»00 


It  should  be  noted  that  equations  8  and  9  contain  no  mention  of  the  enemy's 
initial  commitment  (x  )  or  of  the  rate  at  which  additional  submarines  are 
committed  (r). 

DISCUSSION 

This  model  involves  aggregated  parameters ,  each  of  which  has  an  obvious 
physical  interpretation.  The  appropriate  values  to  use  in  any  particular  case 
must  be  determined  by  separate  analysis. 

The  encounter  rate  (b)  depends  on  whether  the  submarines  are  armed  with 
missiles  and  torpedoes  or  torpedoes  only.  For  a  missile  submarine  "encounter" 
occurs  whenthe  carrier  has  been  localized  closely  enough  that  the  commanding 
officer  feels  he  has  enough  information  to  warrant  firing.  For  a  torpedo  submarine, 
encounter  does  not  occur  until  the  submarine  has  approached  the  screen  closely 
enough  to  risk  detection.  3 


If  there  is  no  screen,  "encounter"  occurs  when  the  submarine  reaches  firing 
position. 
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The  encounter  rate  also  depends  on  the  amount  of  information  available 
to  the  submarines.  If  they  must  search  at  random  for  the  carrier,  an 
encounter  rate  can  be  derived  from  the  average  detection  range,  the  speeds 
of  carrier  and  submarines,  size  of  the  operating  area,  and  the  weapon  firing 
range.  The  detection  range  depends  on  how  much  electromagnetic  and 
acoustic  noise  is  radiated  by  the  carrier  and  its  close  escorts.  If  the  carrier 
if  very  quiet  and  can  be  identified  only  by  visual  sighting,  it  may  be  appropriate 
to  take  account  of  data  suggesting  that  submarines  have  a  near-zero  capability 
in  this  regard  during  darkness. * * 5  6  (It  may  also  be  better  for  missile  submarines 
to  fire  torpedoes  in  this  case. )  If  the  submarines  have  some  external  intel¬ 
ligence  so  that  their  search  is  not  random,  the  encounter  rate  must  be  selected 
to  reflect  the  details  of  the  particular  situation.  For  practical  purposes,  if 
the  carrier  or  its  close  escorts  are  very  noisy  so  that  the  submarines  can 
achieve  very  long  average  detection  ranges  with  their  own  sensors,  it  is  the 
same  as  if  they  had  external  intelligence.  The  encounter  rate  in  such  cases 
may  be  limited  by  a  submarine's  ability  to  maneuver  and  its  boldness,  rather 
than  by  its  ability  to  localize  the  carrier. 

Attrition  rate  a  results  from  the  action  of  all  the  ASW  forces  used  in  any 
particular  case.  It  is  die  rate  of  kill  that  occurs  in  the  absence  of  successful 
attacks  on  the  carrier;6  hence,  it  incorporates  the  kills  that  occur  when  sub¬ 
marines  encounter  the  carrier  and  try  unsuccessfully  to  attack  it,  but  excludes 
kills  resulting  from  successful  attacks.  7  It  can  be  seen  that  the  attrition  rate 
depends  on  the  encounter  rate,  since  each  encounter  offers  some  risk  to 
the  submarine  as  well  as  the  carrier.  The  remaining  component  of  a  is  the 
rate  (a)  at  which  surviving  submarines  are  attrited  in  the  absence  of  encounters 
with  the  carrier.  The  attrition  rate,  then,  also  includes  kills  by  ASW  aircraft 
and  other  systems  that  may  be  present  in  a  particular  case. 

The  attrition  rate  is  one  of  three  parameters  in  this  model  that  incorporate 
estimates  of  ASW  force  effectiveness.  The  others  are  the  encounter  rate 
(in  their  encounter  attempts  submarines  may  be  slowed  down  and  otherwise 
inhibited  by  effective  ASW  forces),  and  the  probability  of  successful  attack 
given  an  encounter  (h). 

4~  ~ 

Reference  (a),  page  7  ff. 

5  References  (b)  and  (c). 

This  definition  of  a  is  a  consequence  of  the  way  equation  1  is  derived 
(see  appendix  A). 

7 

Kills  resulting  from  successful  attacks  are  accounted  for  in  equations 

6  through  9. 
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Values  for  h  can  be  obtained  from  analysis  of  a  submarine's  attack  after 
it  has  encountered  the  carrier.  For  a  missile  submarine,  h  involves  AAW 
as  well  as  ASW  considerations.  For  a  torpedo  submarine,  h  involves  screen 
penetration  probabilities  and  carrier  evasion  probabilities.  The  probability 
that  the  submarine  is  killed  given  that  it  makes  a  successful  attack  also 
involves  different  factors  for  the  two  different  types  of  submarines. 

A  convenient  way  to  summarize  the  possible  outcomes  of  an  encounter, 

no  matter  which  kind  it  is,  is  as  follows:  ^ 

Carrier  hit 


Submarine  killed 


Yes  No 


Yes 

No 


P1 

P2 

P3 

P4 

In  every  encounter,  the  carrier  is  either  hit  or  not  hit  and  the  submarine  is 
either  killed  or  not  killed.  8  Thus: 


Pj  =  Pr  (carrier  hit  and  submarine  killed) 

P2  =  Pr  (carrier  not  hit  and  submarine  killed) 

P3  =  Pr  (carrier  hit  and  submarine  not  killed) 

P4  =  Pr  (carrier  not  hit  and  submarine  not  killed). 

P1  +P2+P3+P4  =  L 


Parameters  h,  a,  and  S,  required  by  the  model,  are  related  to  these 
probabilities  as  follows: 


a  =  bP^  +  a 
S  =  b(P1  +  P2)  +  a. 

The  number  of  successful  attacks  (N)  a  carrier  can  absorb  before  being 
forced  to  retire  depends  on  the  number  of  hits  per  successful  attack  and  the  number 
of  hits  the  carrier  can  take.  These  numbers  are,  in  general,  not  the  same  for 
missiles  as  for  torpedoes.  In  any  event,  N  will  probably  be  based  as  much 
on  judgment  as  on  analysis.  N  can  also  be  varied  to  study  the  question  "what 
would  happen  if  the  carrier  were  constructed  so  as  to  be  able  to  absorb  1,  2  etc 
successful  attacks.  " 

o  — - — 

These  kills  include  all  prosecutions  triggered  by  the  encounter. 

9 

References  (d),  (e),  (f),  and  (g)  contain  data  on  the  vulnerability  of  carriers 
to  missile  or  torpedo  hits. 
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In  summary,  the  parameters  of  this  model  are  determined  by  — 

The  number,  type,  and  effectiveness  of  the  enemy's  submarines, 
The  intelligence  available  to  the  enemy's  submarines, 

The  number,  type  and  effectiveness  of  the  ASW  forces, 

The  details  of  the  situation  and  the  tactics  employed,  and 
The  vulnerability  of  the  carrier. 

As  these  change,  they  cause  the  parameters  of  the  model  to  change.  There¬ 
fore,  the  model  can  be  used  to  study  the  ways  in  which  the  type  of  enemy  sub¬ 
marine  etc.  affect  carrier  on-station  probability  and  enemy  submarine  losses. 

Appendix  B  gives  some  numerical  examples  to  illustrate  the  application  of 
this  model. 
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APPENDIX  A 


APPENDIX  A 


MATHEMATICAL  DERIVATIONS 


This  model  treats  the  carrier-submarine  problem  by  using  differential 
equations.  Five  constant  quantities  described  as  "rates''  are  used  in  formulating 
the  equations:  a,  b,  a,  S,  and  r.  In  some  instances,  it  is  useful  to  think 
these  rates  as  "probabilities  per  unit  time"  that  certain  events  will  occur. 

In  particular,  when  deriving  numericalvalues  of  the  parameters  from  data  or 
theoretical  analyses,  it  is  often  convenient  to  think  in  terms  of  probabilities  and 
expected  values,  rather  than  rates.  The  conversion  of  probabilities  to  rates  and 
vice  versa  is  a  satisfactory  procedure,  provided  one  selects  a  unit  of  time 
small  enough  that  there  is  little  chance  of  an  event's  occurring  more  than  once 
during  a  unit  of  time,  and  yet  large  enough  that  there  is  some  probability  of  its 
occurring  once. 


Let  x(t)  be  the  expected  number  of  attacker's  submarines  alive  at  t,  given 
the  carrier  has  not  been  successfully  attacked  up  to  time  t.  (The  other  symbols 
used  in  this  appendix  are  defined  in  the  main  body. ) 

If  the  carrier  has  not  been  successfully  attacked  up  to  time  t,  the  probability 
that  it  will  not  be  successfully  attacked  during  an  additional  increment  of  time 
dt  is  given  by 


PQ  (dt  1 1)  =  1  -  b  x(t)  h  dt, 


(A-l) 


so  that 

and 

This  leads  to 


P  (t  +dt)  =  P  (t)P  (dt  1 1) 


dPo(t) 


PQ(t)  [1  -  b  x(t)  hdt  ] 


=  P  (tixit)  -  P  (t)  =  -bhx(t)  P  (t)  dt 


P  (t) 

o 


f 


dP 


=  -bh 


; 


x(t)dt 


and 


U 

PQ(t)  =  exp  I  -bh  J  x(t)dt 

(  o 


(A-2) 


(A -3) 


where  f*"  x(t)  dt  is  simply  the  total  number  of  search  days  achieved  by  all 

■'o 

submarines  up  to  time  t,  given  the  carrier  has  not  been  successfully  attacked. 
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The  expression  for  the  expected  number  of  submarines  present,  x(t),  may 
be  written  as 


dx  _ 

at  -  ‘ax  +  r  ■ 


so  that  x(t)  =  r/a  +  e~at  (x  -  r/a)  . 

o 

When  this  expression  for  x(t)  is  used  with  equation  A -3,  the  result  is 

-at\ 


d  /t-\  -bh 

Po(t)  =  exP  ~ 


[rt  +  (xQ  -  r/a)  (l-e“ac)  ]  . 


(A -4) 

(A  -5) 

(A -6) 


An  important  special  case  occurs  when  r=0,  meaning  that  the  enemy 
commits  no  additional  submarines  after  the  action  begins.  In  this  case 


PQ(t)  =  exp 


t- 


bhx 


o  -at- 

-  (1-e  ) 


(A -7) 


When  r=0  it  is  possible  that  all  xq  of  the  submarines  may  be  killed  before 
the  carrier  is  successfully  attacked.  The  probability  of  this  event  is 

bhx 


lim  PQ(t)  =  exp 
t  “ 


(A-8) 


An  expression  for  P  (t),  the  probability  that  the  carrier  has  been  successfully 

attacked  exactly  n  times  by  time  t,  is  derived  as  follows.  Denote  by  P  (t-T) 

the  probability  that  there  is  no  successful  attack  in  the  interval  (T,  t),  given 
that  the  carrier  is  on  station  at  T.  Then  P  (t-T)  is  given  by  equation  A-6 

except  that  the  xq  of  equation  A -6  is  replaced  by  x(T),  and  t,  by  t-T;  that  is, 

relative  to  the  interval  (T,  t),  the  initial  number  of  submarines  is  x(T). 


Intuitively,  it  would  be  expected  that 


Pq  (t-T)  Pq  (T)  =  Pq  (t)  , 


(A-9) 


and  by  use  of  equation  A-6  it  is  easy  to  show  that  this  is  correct.  Now  for 
n  <  N  , 


5n<')=f 


f 2 

r  -| 

r  (A-10) 

tj=0  Po<'-'„> 

bhx(t  )  dt  P(t  -t  .) 
n7  n  n  n-l/J 

•  •  • 

bhxft^dtj  Po(tj) 
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Using  equation  A-9  reduces  this  to 


jyt)  ~  p0(t)  5  '•••  5  2  bhx(t  )  dt  .  ..bhx(t  )  dt, 

t  =0  t,=0  n  n 

n  1 

From  equation  A -2 

i  dP  (t  )  d(ln  P  (t.) ) 

bhx(t.)=  -*±7-,  ■  *  =  -  _ 0  1 

1  W  dti  dt. 

Using  equation  A-12  in  equation  A-ll  and  integrating  gives 

n 


pn(t)  =  pQ(t) 


t-ln  PQ(t)  ] 
W~ 


(A-ll) 


(A-12) 


(A-13) 


Since  this  is  the  probability  of  exactly  n  successful  attacks,  the  probability 
of  N  or  fewer  successful  attacks  by  time  t  is  given  by 


TT 


N 

E 

n=0 


p  (t)  . 
n 


(A -14) 


Next  an  expression  is  derived  for  K(t),  the  expected  number  of  sub¬ 
marines  killed  up  to  time  t.  From  the  definition  of  the  process  that  is 
modeled  here, 


dK(t) 

dt 


rrN(t)  S  x(t)  . 


(A-15) 


By  use  of  equations  A-12,  A-13,  and  A-14  this  can  be  written  as 


dK 

dt 


S_ 

bh 


dP 
_ o 

dt 


N 

E 

n=0 


L-ln  P  I1 
_ _ o  J 

n ! 


Therefore,  K(t) 


N 


p0(t) 

s 


^-ln  Po-)n 

n! 


dP 


o 


(A -16) 


(A-17) 
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Carrying  out  the  integration  of  an  individual  term  gives 


Pol') 


[-In  Pi" 


n  ! 


dP  =  TT  (t) 
o  n  w 


so  that  K(t)  = 


^  n=0  U  '  J  ' 


(A-18) 


(A-19) 


The  "exchange  ratio" 


R(t) 


K(t) 


1- 


TT 


N 


R(t) 
(t)  • 


is  defined  by 


For  N=0, 

R(t) 


_s_ 

bh 


for  all  values  of  t .  For  N  >  0  , 


(A-20) 


(A -21) 


lim  R(t)  =  (N+l)  ,  (A -22) 

t  -*00 

provided  that  r  >  0  .  The  interpretation  of  equation  A-22  was  given  in  the 
main  body. 


ALTERNATE  WAYS  OF  VIEWING  THE  CARRIER -SUBMARINE  INTERACTION 
PROBLEM 

Discrete  Time  Intervals 


The  carrier-submarine  interaction  problem  has  been  frequently  modeled 
by  considering  fixed  intervals  of  time  At  so  that  the  "time  variable"  became 
the  number  of  time  intervals  (t/At) .  The  various  appropriate  probabilities 

were  then  raised  to  this  power:  such  as  (1  -  aAt)^^  .  It  may  be  shown  that  the 
analytical  results  of  such  a  model  transform  directly  to  the  results  derived 
earlier  in  this  paper  as  the  time  interval  At  is  allowed  to  become  arbitrarily 
small. 
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Unconditional  Attrition  Rate  Constant 


It  would  seem  that  a  model  could  be  developed  which  assumed  the  un¬ 
conditional  attrition  rate  of  submarines  as  constant.  This  has  been  done.  1 
The  results  are  the  same  when  the  carrier  is  not  very  likely  to  be  hit  and  some¬ 
what  more  optimistic  than  predicted  in  the  body  of  this  paper  when  the  carrier 
is  likely  to  be  hit.  However,  this  development  requires  the  simultaneous 
solution  of  many  differential  equations  and  the  solutions  for  any  but  the 
simplest  cases  are  complex  and  difficult  to  evaluate.  A  generalization  of  this 
development  to  the  case  of  many  carriers  operating  in  the  same  area  has  also 
been  completed.  2 


J.  K.  Tyson,  appendix  C. 

2 

Reference  (h) . 
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SAMPLE  NUMERICAL  CALCULATIONS 

This  appendix  illustrates  how  tactical  data  can  be  used  to  get  numerical 
values  for  the  parameters  of  the  model.  Sample  calculations  are  made,  and  the 
resulting  on -station  probabilities  and  exchange  ratios  are  displayed.  To  keep 
the  presentation  unclassified,  the  numbers  used  are  arbitrary;  they  were  not 
derived  from  exercise  data  or  theoretical  studies.  1 

Cases  of  an  escorted  carrier  opposed  by  torpedo  submarines  and  by 
missile  submarines,  and  an  unescorted  carrier  2  opposed  by  torpedo  sub¬ 
marines  are  illustrated.  An  escorted  carrier  can  be  detected  at  longer  ranges, 
owing  to  the  active  sonar  of  its  escorts;  encounter  rates  reflect  this  increased 
detection  range.  An  unescorted  carrier  can  be  attacked  immediately  once  it 
has  been  encountered,  with  virtual  assurance  that  the  attack  will  be  successful. 

The  arbitrary  probabilities  used  to  describe  the  four  possible  outcomes 
of  submarine-carrier  encounters  are  given  in  table  B-I. 


TABLE  B-I 


C  arrier 
hit 

Submarine 

killed 

Yes 

No 

Escorted  carrier  vs. 

torpedo  submarines 

Yes 

.21 

.29 

No 

.09 

.  41 

Unescorted  carrier  vs. 

torpedo  submarines 

Yes 

.  15 

0 

No 

.  85 

0 

Escorted  carrier  vs. 

missile  submarines 

Yes 

.05 

.  10 

No 

.28 

.  57 

Encounter  rates  were  calculated  by  methods  described  in  reference  (a),  an 
unclassified  document. 

2 

Operating  at  slow  speed  and  under  strict  electronic  silence  in  a  dispersed 
formation. 


Escorted  carrier  versus  torpedo  submarines: 

b  =  encounter  rate  per  day 

=  •  25  (corresponds  to  a  detection  range  of  35  miles  in 
an  operating  area  of  70,000  dquare  miles) 

a  =  b  •  Pr  (submarine  killed  given  that  the  carrier  is  not 
hit)  +  rate  of  kill  from  all  other  causes 

=  .25  •  .  29  +  .  05  =  .  12 

h  =  Pr  (carrier  hit,  given  an  encounter) 

=  .  21  +  .  09  =  .3 

S  =  unconditional  rate  of  submarine  attrition 
=  .25  (.21  +  .29)  +  .05  =  .18  . 

Unescorted  carrier  versus  torpedo  submarines: 

b  =  .  02  (corresponds  to  a  detection  range  of  7  miles  in  an 
area  of  70,000  square  miles)3 

a  =  .  02  (0)  +  .  05  =  .05 

h  =  1.0 

S  =  .  05  . 


Unescorted  carrier  versus  cruise  missile  submarines:  If  the  carrier  is 
unescorted  the  submarine  must  see  it  to  find  it  and  should  thus  torpedo  it.  This 
is  therefore  identical  to  the  previous  case. 

Escorted  carrier  versus  cruise  missile  submarines: 

b  =  .25 

a  =  .  25  •  .  10  +  .  05  =  .  075 

h  =  .33 
S  =  .09  . 


Carrier  average  SOA  10  knots,  submarine  average  SOA  3  knots, 
marine  can  identify  the  carrier  only  during  daylight. 


The  sub- 
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The  numerical  results  of  these  three  examples  are  given  in  table  B-II, 
for  xq  =  5  and  r  =  0  in  each  case.  The  probability  of  zero  successful 

attacks  up  to  time  t,  PQ(t),  is  shown  in  figure  B-l.  The  probability  of 

the  carrier’s  being  on  station  at  time  t,  tTj  (t),  is  shown  in  figure  B-2,  for 

hi  1  a 

TABLE  B-II 

NUMERICAL  RESULTS  OF  SAMPLE  CALCULATIONS 


Escorted 

carrier 

versus 

torpedo 

submarines 

Unescorted 

carrier* 

versus 

torpedo 

submarines 

Escorted 

carrier 

versus 
missile 
s  ubmarines 

Encounter  rate,  b 

.25 

.  02 

.25 

Conditional  submarine  attrition  rate,  a 

.  12 

.  05 

.  075 

Overall  submarine  attrition  rate,  S 

.  18 

.  05 

.  09 

Pr  (successful  attack  given  encounter),  h 

.  3 

1.  0 

.  33 

Exchange  ratio  (if  N=0),  R 

2.  3 

2.7 

1. 1 

Exchange  ratio  (if  N=l),  R 

4.7 

5.  3 

2. 1 

Pr  (carrier  on  station  after  10  days, 

N=0),  tto(10) 

.  12 

.  45 

.06 

Pr  (carrier  on  station  after  10  days, 

N=l),  tt1(10) 

.  36 

.  81 

.21 

*  "Quiet"  carrier  in  dispersed  formation. 
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FIG.  B-1:  ON-STATION  PROBABILITY  OF  A  CARRIER  THAT  MUST  RETIRE  THE 
FIRST  TIME  IT  IS  SUCCESSFULLY  ATTACKED 


FIG.  B-2:  ON-STATION  PROBABILITY  OF  A  CARRIER  THAT  CAN  ABSORB  ONE 
SUCCESSFUL  ATTACK  AND  STILL  REMAIN  OPERATIONAL 
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A  MODEL  OF  CARRIER -SUBMARINE  INTERACTIONS 
WITH  UNCONDITIONAL  ATTRITION  RATE  CONSTANT 

INTRODUCTION 

The  body  of  this  paper  describes  a  comparatively  simple  model  for 
carrier-submarine  interactions  which  provides,  as  a  function  of  time,  a 
probability  distribution  for  the  number  of  successful  attacks  on  the  carrier 
and  the  expected  number  of  submarines  surviving  the  engagement.  This 
appendix  addresses  the  same  basic  model,  but  follows  a  different  procedure 
which  avoids  certain  difficulties. 

In  brief,  the  assumptions  of  this  development  are  as  follows: 

a.  Each  submarine  in  the  area  will  independently  and  randomly 
encounter  the  carrier  and/or  its  local  defenses  at  a  constant  rate  denoted  bv 
"b". 


b.  Each  encounter  will  result  in  one  of  four  possible  events,  with 
associated  probabilities: 

1.  With  probability  P  ,  the  carrier  is  hit  and  the  submarine  is 

eliminated.  i 

2.  With  probability  P2,  the  carrier  is  not  hit  and  the  submarine  is 

eliminated.  z 

3.  With  probability  P„,  the  carrier  is  hit  and  the  submarine 
survives  the  encounter. 

4.  With  probability  P.,  the  carrier  is  not  hit  and  the  submarine 

survives  the  encounter.  ’ 

c.  Each  submarine  will  independently  and  randomly  encounter  and  be 
eliminated  by  area  defenses  at  a  constant  rate  denoted  by  "a". 

d.  Additional  submarines  may  enter  into  the  engagement  at  random  with 
a  constant  rate  denoted  by  ”r". 

These  assumptions  define  a  continuous  time  Markov  chain  in  which  the  state 
of  the  system  is  described  by  the  number  of  "hits"  received  by  the  carrier  and 
the  number  of  surviving  submarines  in  the  engagement.  The  state  probabilities 
are  denoted  by  ffl  h(t)  ,  defined  as  the  probability  at  time  t  that  the  carrier  has 

received  h  hits  and  is  engaged  by  n  surviving  submarines.  The  possible 
transitions  between  states,  and  their  associated  rates,  are  shown  in  the  following 
table:  s 


C-l 


Transition  rate 


State  change 


n,  h 

n-1,  h+1 

nbP: 

n,  h 

n-1,  h 

n  (b  P 

n,  h 

n,  h  +  1 

n  b  P 

n,  h 

n  +  1,  h 

r 

The  state  probabilities  are  determined  by  solving  the  following  set  of 
differential  equations: 


df 


n,  h 
dt 


=  -n  [«P1+P2+P3)  +  or]  fn  h  -r  fnh 

■Kn+DbPj  fn+iih.j  +  <n+l)[bP2  +  a]  fn+1_ 


+  nbPQf  +  rf  .  , 

3  n,  h-1  n-l,h  . 


(C-l) 


This  equation  applies  for  all  n,  h,  with  the  understanding  that  terms  involving 
f  with  negative  subscripts  are  to  be  suppressed.  The  initial  conditions,  when 
N  submarines  are  engaged  at  the  start,  are 


fn  k  -  0  all  other  n,  h  . 

SPECIAL  CASE,  ONE  OPPOSING  SUBMARINE 


The  set  of  equations  C-l  is  simplified  somewhat  if  only  one  submarine 
is  initially  present  and  there  is  no  reinforcement.  Then  only  the  state 
probabilities  f^h  and  fQ  h  need  be  considered,  and  the  set  C-l  specializes 


df 


1,0 


dt 


df 


l,h 


dt 


-[b(Pi  +  P2  +  P3)  +  or]  fl  o 

-lb(Pl  +  P,  +  P3)  +  o]  f1>h+W>3tlih.1 


(C-2. 1) 


C-2 


(C-2.2) 


df 


^  =  <bP2  +  “>  fl>0 


df 


dt  bPl  fl,  h-1  +  ^bP2  +  ^  fl,  h 


The  initial  conditions  are 


f,  =  1 
1,  o 


fl,h  fo,o  fo,h  =  0 


(C-3) 


The  probability  that  the  carrier  receives  no  hits  is,  by  definition,  the 


sum  of  f.  and 

1,0 

C-2.  1  and  C-2.  2: 

f  .  To  determine 

0,  0 

fl,o 

and  fQ  Q  we  integrate  equations 

fl,  0 

exp  -  |_b  (P1  +  P2 

+  p3> 

+  a  ]  t 

f  - 

bP2  +  a  1 

1-exp 

-  [b(P1  +P2  +P3)  +  a~]  t^  . 

0,  0 

b(Pl+P2+P3>  +  a  V 

The  probability  the  carrier  receives  no  hits  is  now  given  by 
b(P1+P3) 

Po  =  1  "  b(P1+P2+P3)  +a 

If  N  submarines  are  initially  present  and  are  not  reinforced,  we  would  expect 
the  probability  the  carrier  receives  no  hits  to  be  given  by 


0 


exp 


-  b  (P^  +  P2  +P3)  + &  ]  t 


) 


P 

o 


UP1+P3) 

j1  “  b(p1  +P2  +  P3)  +c* 


^1-exp-  [b(P1+P2+P3)  +a]t 


)i 


N 


in  view  of  the  independence  assumption  of  the  model.  This  result  differs  from 
equation  2  of  the  main  body  in  functional  form  presumably  because  the  latter  is 
based  at  the  outset  on  an  expected-value  calculation  of  the  number  of  submarines 
versus  time.  Of  possibly  greater  significance,  the  parameter  combination 

h^Pl  +P2  +IV  +  a  here  reP*aces  the  expression  a  =  bP2  +  a  used  in  the  main  body. 
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To  find  the  distribution  of  the  number  of  hits  received  by  the  carrier 
requires  the  recursive  integration  of  the  full  set  of  equations  C-2.  This  is 
painlessly  accomplished  by  introducing  the  probability  generating  function,  1 
defined  as  follows: 


VX.t>  =  E  f i _ h  <t)  Xh 


F0<X’t)  =  S  xh  • 
h=0  ’ 

From  equations  C-2  and  C -3  it  follows  that 


-  [b(Pj  +P2  +P3)  +  <*]  F2  +  X  bP3F1 


(bP2  +a)  F1  +  X  bP1F1 


d^ 

dt 


dF, 


dt 


Fj  (X,0)  =  1 
Fq  (X,0)  =  0  . 

The  solution  of  these  equations  is  easily  seen  to  be 

Fl  t)  =  exP  {-  [b(P1+P2+P3)  +  a]  t  +  X  bP3t  } 


bP9  +a  +  XbP. 

F0  (X’  C)  ~  b(Pi+P2+P3)+«-Xbp3 


[l  _exP  l  “  [b(pi  +p2 p3)  +a]  +  XbPgt  }] 


jHie  generating  function  for  the  probability  that  the  carrier  receives  exactlv 
hits  is  just  the  sum  of  Fj  and  FQ.  If  N  submarines  are  present  and  are 

not  reinforced,  we  would  expect  the  generating  function  to  be  (F3+F0)N  ,  in  view 

of  the  independence  assumption  of  the  model.  Further,  the  generating  function 


1  jy*,  F?ller’  ”An  Intr°duction  to  Probability  Theory  and  its  Application, 
Vol.  1,  p.  264. 
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for  the  probability  that  the  carrier  receives  h  or  fewer  hits  is  given  by^ 


"  LFpX.t)  +F0(X,t)jN 

S  ^h  (t)x  =  - Px - 

0 


where  the  left-hand  side  employs  the  notation  given  in  the  main  body.  This 
result  reduces  the  problem  of  finding  the  tt^'s  to  a  routine  but  tedious  process 

of  power  series  expansions.  For  example,  when  N  =  1,  and  in  the  limit  t  -*  ®  , 

bP^  +  a 

P0=n0  =  b(P1+P2+P3)  +a  ’ 

and  if  h  >  1  , 

[b(P1+P2)  +  a]  [b(P1+P3)  ]  (bP3)h_1 
h  [b(P1+P2+P3)  +a]  h+1 


1  "^b(P1+P2+P3)  +  a  j 

Finally,  the  probability  that  the  submarine  survives  for  time  t  is  given 
by  Fj  (1,  t)  =  exp  -  [b(P1+P2)  +  a  ]  t  . 

GENERAL  CASE 


The  solution  of  the  more  general  equations  C-l  can  be  carried  out  by  a 
process  similar  to  that  described  above.  A  generating  function  for  the  joint 
distribution  of  the  two  random  variables  n  and  h  is  defined  as  follows: 

CO 

F(Z,  X,  t)  =  S  f  7nXh  . 

n,  h 

n,  h=0 

The  equation  governing  F  is  determined  by  multiplying  each  side  of  equation 


2 


Op.  cit. ,  p.  265. 
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C-l  by  ZnXh  and  summing, 
as 


The  final  result,  on  making  use  of  identities  such 


nZ 


n  _ 


Zaz11 

’ 


is  found  to  be 


|f  +  Kl-ZJF 


+  l  [(P1+P2+P3(l-X))  +aj  Z-  ^(PjX+P^  +aj  }  ||  =  o  , 

(C-4) 

with  initial  condition,  describing  the  presence  of  N  submarines  and  no  hits 
received  by  the  carrier,  given  by 

F(Z, X,  0)  =  ZN  .  (c-5) 

For  convenience,  the  following  "dimensionless"  parameters  are 
introduced: 

t  =  [b(P1+P2+P3)  +a-  XbP3  jt 

P  =  TMPprp-Fpp-^rxEFyy  (c.6) 


bP2  +  a  +  XbP1 
CT  ~  b  (P1+P2P3)  +  -  XbP3 

Equation  C-4  then  takes  the  form 

§f  +  p  (l-Z)F  +  (Z  -  „)  g  =  0. 

A  "steady  state"  solution  is  defined  by 
p(l-Z)  Y+(Z  -a)  &  =  o 

or  Y(Z)  =  exp  -  f  p  d  Z’ 

J  v  L  -  (j 

=  (Z  -  a)~  P  ep  (Z-cr) 
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Equation  C-7  is  transformed  to  a  more  convenient  form  by  making  a 
substitution,  F(Z,t)  =  Y  (Z)  G  (Z,t)  ,  which  leads  to  the  following  equation 
for  G: 

+  (Z  -  a)  H  =  0  .  (08) 

The  general  solution  of  equation  C-8  is  an  arbitrary  function  of 


or,  equivalently,  a  function  of  (Z-g)e  T.  Thus,  the  general  solution  of 
equation  C  -  7  is  given  by 

F  (Z,  T)  =  (Z  -  CT)"P  e"P  (Z_cr)  Gt  e  T  ]  . 


The  functional  form  of  G  is  established  by  making  use  of  the  initial  condition, 
equation  C-5.  The  result  of  this  process  is  the  solution 


F  (Z, X, t)  =  e  pr(Z-rr)(l-e  T)  -  (1-CT)  t  ]  [ CT+(z_CT)e-r  ]N  ^ 


(C-9) 


vdiere  the  parameters  p,  <j,  T  are  functions  of  X  as  defined  by  equations 

C-6. 

Particular  results  can  be  derived  from  C-9  as  desired.  For  example,  the 
expected  number  of  submarines  at  time  t  is  given  by 

n  =  ^F  I  X=1 

BZ  I  Z=1 


_ -  1  _P  "  [  b(Pi  "^P9 )  +^1 1 1  ,  xt  -  (" t>(p i  +p0 )  +Q'l  t  . 

b(P~fTpT«  |_  1  2  J+Ne  L  1  2  J 

The  distribution  of  the  number  of  hits  received  by  the  carrier  has  the  generating 
function  given  by 


F(l,  X,  t)  = 


bp2+a+XbP 


1  +b(p1+p3)  ( 1  -  X)  e  “[b^Pi+P2+P3)+Q']t  +  XbP3t'|N 
’  +po)  +a  -  XbP0  J 


MPyPp^+p^y  ,  - 

2 


rb  (P^  +p3)  (1-X) 

X6XP  rblPj+^T^)  +a-XbP3) 


1  "[b(p1+p2+p3)+^]t  +  Xbp3t-e  -[b(Pi+P2+p3)+^]t+Xbp3t 


(C-10) 
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The  probability  that  the  carrier  receives  exactly  h  hits  is  obtained  by 

finding  the  coefficient  of  X  in  the  right-hand  side  of  equation  C- 10.  The 
calculations  are  straightforward  but  lengthy.  For  convenience  in  presentation 
of  results,  the  probability  distribution  will  be  given  as  the  convolution  of  two 
distributions : 

1.  The  probability  of  h  hits,  given  that  N  submarines  are 
engaged  initially  and  there  is  no  reinforcement:  Ph  (t  |  N,  0)  . 

2.  The  probability  of  h  hits ,  given  that  no  submarines  are 
engaged  initially  and  that  submarines  are  introduced  at  a  rate  r  : 
Ph  (t  1 0,  r)  . 

The  following  definitions  are  employed: 


X  = 


P  ~ 


t  = 


Q,(t)  = 

k=o 


b(P!+P3> 

+p2  +P3)  +  a 

rb  (Pj  +P3) 

Lb<P1+P2+P3)+aJ2 

bP3 

b(P1+P2+P3)  +a 
[b  (p1+p2+p3)  +a]  t 


The  quantities  Qj  are  the  "tails"  of  a  Poisson  distribution  with  a  parameter 
r  ;  equivalently,  they  are  the  probability  that  the  sum  of  (J+  1)  exponentially 
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distributed  variables  has  a  value  less  than  t  .  With  these  definitions, 
for  the  first  few  P^'s  have  the  following  explicit  forms: 

P0(t|N,  0)  =  (1-XQq)N 

Px  (t  IN,  0)  -  nX(i-\q0)n-1(q0-Cq1) 

P2(t|N,  0)  =(2)^2(1^Q0)N'2(Q0-CQ1)2+NXC(1-XQ0)N'1(Q1 
p3(t|N,  0)  =  (1-A.Q0)N-3  (Qq-Cq/ 

+  2  (2  )x2  c  (1_xqo)N"2  (QrCQ2) 

+  nxc2(i-xq0)n_1  (Q2-Cq3) 

p4  (t  I  N,  0)  =(4)^4(1^Q0)N_4+3(3)x2  C(l-XQ0)N_3(Q0-CQir 
+  (^)x2C2(1-XQ0)N‘2{(Qi-CQ2)2  +2(Q0-Q1)(Q2- 

+  NXC3  (I-^Qq)^1  (Q3-CQ4) 


PQ  (t  |  0,  r) 
Pj  (t  10,  r) 

P2(tl0,  r) 

P3  (t  10,  r) 


e-P(-Q0) 

e"P(T"Q0)  p[(l-C)TQ0  -(1-2  0QJ 


,-P  (T"Qr 


2  f 
o){fr  (1'°  tQq-o-zOq^ 


+  pC 


(I-C)tQ 


,-P  (T“Qn)J  P 


0*1^  (1-0  TQ0  -(1-20  qJ" 

+  P2Z  [(l-o  tQq-(1-20Q1]  [(1-OtQ1  -  (2-3C)  Q2] 
+  pC2  [(l-o  TQ2  -  (3-40  Q3]| 


the  results 

-cq2) 


■(Qj-CQj) 

Vj 


j-(2-30Q2J| 
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In  view  of  the  complexity  of  these  expressions,  it  seems  worth  pointing 
out  one  special  case.  If  it  is  justified  to  assume  that  the  carrier  can  survive 

h  hits  with  probability  q  ,  then  a  comparatively  simple  expression  for  the 
probability  that  the  carrier  survives  for  time  t  is  provided  by  equation  C-10 
with  the  input  X-q.  In  the  notation  used  directly  above  the  result  is 


ho-cq)  -  i+e"T(1‘?<i) 


(l-Cq) 


Probability  of  survival  =  e 


This  result,  of  course,  can  be  derived  directly  from  the  basic  model  with 
appropriate  modification  of  the  probabilities  describing  the  result  of  an 
engagement. 

POISSON  APPROXIMATION 

Returning  to  equation  C-10  and  the  dimensionless  parameters  defined  after, 
we  note  that  if  the  parameter  C  is  set  equal  to  zero  the  probability  distribution 
reduces  exactly  to  a  convolution  of  Poisson  and  binomial  random  variables. 

This  parameter  appears  to  account  for  multiple  successful  attacks  by  an  individual 
submarine.  This  observation,  together  with  the  fact  that  a  binomial  distribution 
can  be  approximated  by  a  Poisson  distribution,  suggests  that  the  total  hit 
distribution  can  be  approximated  by  a  Poisson  distribution  -  indeed,  this  is  the 
result  derived  in  this  paper.  There  are  various  ways  to  obtain  such  a  rep¬ 
resentation;  the  most  straightforward  is  to  use  the  generating  function,  equation 

The  generating  function  is  written  in  the  form  F(l,  X,t)  =  exp  3>(X),  where,  of 
course,  $(X)  is  the  logarithm  of  the  right-hand  side  of  equation  C-10.  The 
coefficients  are  defined  by  the  power  series  expansion  of  $>  (X).  Thus, 


C-10 


F(l,X,t)  =  exp  jfl>0  +  3^  X  +  3>2  X2  +  . . .  ]  . 


The  next  step  is  to  expand 
after  some  rearrangement  the  result  becomes 


exp  $2  X2  +  3>3  X3  +  . .  .J  as  a  power  series; 


F(  1 ,  X, t)  =  e$0  +  $1  e$l  (X_1) 


3> 


1  +3>2  X  2+  3^  X  3+  ( — -  +3>4)X4+... 

(C-ll) 


Equation  C-ll  demonstrates  that  the  distribution  of  hits  can  be  expressed  as 
a  sequence  of  Poisson  probabilities  multiplied  by  a  normalization  factor.  For 
example,  the  cumulative  distribution  of  hits,  tt  ,  can  be  written  as 


TV  =  A 


Sh<4>>  +  *2  Sh-2  <«>>+  *3  Sh-3  <*1'  + 


(C  -12) 


where  A  =  exp  (3>Q+  )  and  are  the  cumulative  probabilities  of  a  Poisson 

distribution  having  a  parameter  .  This  result  is  exact  if  the  expansion  is 

carried  out  to  the  final  term  in  S^;  whether  or  not  the  first  few  terms  are  a 

good  approximation  is  determined  by  the  magnitude  of  the  higher-order 
coefficients  3>2,  3>3  ,  etc.  The  parameters  in  equation  C-12  take  the  following 

form; 


.N  [NX 


A  =  <1-XQ0)  exp-^™-  (Qo-CQj)  -  C  (tQ0  -  2Qp 


3>  = 


NX 


1  1-XQ0  ^0~^CV+P  (T  "  Qq)  "  ^  (TQo  "  2(V 


^2  Q!  "  ^2)  "  /^cL'\2  (Q0_?Q  1)2  +  (TQ0"2Qi)_^TQr3Q2) 


(1  -XQ0) 


V^c2«2-^3,-(-y^2  OQ0-CQ1)(Qi-CQ2)+i^ 


(i-xQ0) 


+ 


(tQ1“3Q2)  -  C(tQ2  -4Q3) 


C-ll 
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From  these  expressions  it  can  be  seen  that  if  NX  <  <  1  and  C  <  <  1,  then 
A  -»  1  and  ^3*  etc.,  -»  0,  and  equation  C-12  reduces  to  the  form  given  in 
the  main  body. 
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